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Three dimensional (3D) Dirac semimetal is a novel state of quantum matter, characterized by the 
gapless bulk four-fold degeneracy near Fermi energy. Soon after its discovery, the classification of 
stable 3D Dirac semimetals with inversion and rotational symmetry have been studied. However, 
only ten out of thirty-two point groups have both inversion and rotational symmetry, and we need 
a more complete classification of stable 3D Dirac semimetals. Here we classify stable 3D Dirac 
semimetals with reflection symmetry and rotational symmetry in the presence of time reversal 
symmetry, which belong to seventeen different point groups. These systems include the systems 
preserving inversion symmetry except C 3 i. They have two classes of reflection symmetry, with the 
mirror plane parallel to rotation axis and the mirror plane perpendicular to rotation axis. In both 
cases two types of Dirac semimetals are determined by four different reflection symmetries. The first 
type of Dirac semimetals will appear through accidental band crossing (ABC). The second type of 
Dirac semimetals have a Dirac point at a time reversal invariant momentum (TBC). We show that 
in both mirror parallel and perpendicular cases, 6 * 2,3 symmetry can only protect stable Dirac points 
via TBC, while 64,6 symmetry can have stable Dirac points as ABC or TBC. We further discuss 
that Weyl line nodes and Dirac semimetal can exist in Brillouin zone at the same time using C 4 v 
symmetry as an example. Finally we classify Dirac line nodes and Weyl line nodes to show in which 
types of mirror plane they can exist. 


I. INTRODUCTION 

Dirac semimetals are new states of quantum mat¬ 
ter. They have gap closing(Dirac points or Dirac line 
nodes) of conduction band and valence band which show 
pseudorelativistic physics of 3-dimensional(3D) Dirac 
fermions near the Fermi energy. Before the discovery of 
Dirac semimetals, the topological quantum states, such 
as topological insulator pj[4] and graphene [5] 6 . can have 
2-dimensional(2D) Dirac fermions. Different from topo¬ 
logical insulators and superconductors, Dirac semimetals 
hold nontrivial features in the bulk states [7Hl4] . When 
the systems reach a quantum critical point between nor¬ 
mal insulator and topological insulator, the accidental 
crossing of inverted bands will generate Dirac points in 
bulk [15] . 

To realize Dirac semimetal we can tune the chemi¬ 
cal composition to the critical point of quantum phase 
transition [15], however those Dirac points are not sta¬ 
ble. Soon stable 3D Dirac semimetals have been the¬ 
oretically predicted [9, 16 and observed experimentally 
in CdsAs 2 and NaqBi |T7ff21] by Angle-resolved photoe¬ 
mission spectroscopy (ARPES). In these materials there 
are two stable Dirac points in k z axis stabilized by ro¬ 
tational symmetry. While in /3-cristobalite structure 
such like BiC^p!], the Dirac points exist at a time re¬ 
versal invariant momentum(TRIM). The unique elec¬ 
tronic band structures of Dirac semimetals make them 
have unusual high mobility, oscillating quantum spin 
hall effect(QSHE) and giant diamagnetism [19]. Re¬ 
cently, a new type of Dirac semimetal, with Dirac line 
nodes(DLN), has been proposed in Cu3NPd|f0j [14] and 


LaN[22|. DLN can exist in the system with or without 
spin-orbit coupling(SOC) [23 . At the same time, theoret¬ 
ical prediction shows that time reversal symmetry (TRS) 
breaking systems including HgCr2Se4[24] has Weyl nodes 
and a Weyl line node(WLN) in its mirror plane. Also sys¬ 
tems with TRS breaking such as pyrochlore iridates[25] 
and (CdO)2(EuO)2|26] or with TRS such as TaAs, NbAs, 
NbP and TaP f27H37] have Weyl nodes. 

Inspired by these work, we ask the question that which 
point group can protect Dirac semimetals. Yang and Na- 
gaosa have classified Dirac points in physical systems pre¬ 
serving inversion symmetry[38]. However, only ten out 
of thirty-two point groups have both inversion and rota¬ 
tional symmetry, and we need a more complete classifi¬ 
cation of stable 3D Dirac points. Meanwhile, considering 
the critical role that reflection symmetry plays in newly 
predicted DLN and WLN semimetals, a classification of 
them by reflection symmetry becomes necessary as well. 
Here we classify 3D stable Dirac points in the systems 
preserving TRS, reflection symmetry and uniaxial rota¬ 
tional symmetry. Known that the reflection symmetry 
plays an important role in classification of topological 
phases[39l Ff0| , we first study the classification of reflec¬ 
tion symmetries in space groups. Apart from point group 
C^i, all the other point groups preserving inversion sym¬ 
metry are covered in our discussion. Then we show that 
four different reflection symmetries can protect two dif¬ 
ferent types of Dirac semimetals, Dirac semimetals cre¬ 
ated by ABC or by TBC. We further discuss the coexis¬ 
tence of Dirac points and Weyl line nodes through k • p 
theory using C 4 V symmetry as an example. Finally we 
classify DLN and WLN to show in which types of mirror 


2 


plane they can exist. 


II. RESULTS 


Consider four energy bands which could be generated 
in a system preserving time reversal symmetry(TRS) 
and uniaxis rotational symmetry, we can describe them 
through a 4 x 4 Hamiltonian in a very general form, 


H = 


h n (k) h n (k) 
h lA) hl$) 



systems, reflection symmetry can be a point group sym¬ 
metry as shown in FIG.[l](ai) or a nonsymmorphic glide 
plane symmetry which is the combination of a reflection 
and a translation t as shown in FIG. l](a 2 ). Neupane et. 
aljTT] realized Dirac semimetal on CdsAs 2 which belongs 
to nonsymmorphic space group IAicd (C^y). Whereas in 
the second class, when k z axis is perpendicular to the 
mirror plane, inversion symmetry will emerge through 
the combination of reflection symmetry and C 2 , C 4 or Cq 
rotational symmetryf38i [41]. Here the point group re¬ 
flection symmetry is shown in FIG. [l](a 3 ) and the glide 
plane symmetry in nonsymmorphic space group is shown 
FIG.0(a 4 ). 
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where cq represents the spin space and t* represents or¬ 
bital space, hfjcj' (a i) is a 2 x 2 matrix and repre¬ 
sent opposite spin in k z direction. All the a^-(fc) are real 
functions and we can determine the parity of each coef¬ 
ficient dij(k) through TRS H(—k) = TiJ(&)T _1 , where 
T = ia y K . The system is invariant under C n rota¬ 
tional symmetry which gives CnH^k^C^ 1 = H(R n k ), 
where R n is a rotation operator for 3D n-fold rotation 
in k-space and the basis are chosen to be eigenstates 
of C n . We choose the rotation axis as k z axis, there¬ 
fore R n k z = k z . Rotational symmetry suggests com¬ 
mutation relation C n H{k z )C ~ 1 = H(k z ) in k z between 
the Hamiltonian H(k z ) and C n operators. Therefore 
we choose a set of bases to make C n a diagonal form 
C n = diag(a p ,aq,a r ,a s ), where a p = exp[i^(p+ §)], p 
is the angular momentum for point group rotational sym¬ 
metry and the effective angular momentum for a screw 
rotation. Here each basis is a C n rotation eigenstate and 
has a definite rotation eigenvalue p+ In the presence of 
TRS, C n has only two independent basis whose rotation 
eigenvalues are a p and a q (see Appendix). 

After applying TRS and rotational symmetry, the 
Hamitonian in k z axis becomes a diagonal form with the 
above basis 

H(k z ) = < 2 oo + &03CT3 + ^ 33 ^ 3^3 + & 3 (U 3 ( 1 ) 

The Dirac points are created only when < 203 , 33 , 30 (^ 2 , m ) = 
0, where m is a control parameter. There are three equa¬ 
tions and two variables, so we need additional symmetry 
constraints to guarantee these equations having at least 
one solution that can generate stable Dirac points. We 
impose reflection symmetry into the systems with uniax¬ 
ial rotational symmetry and TRS, and show the condi¬ 
tions where Dirac points exist. 


II. 1. Classification of reflection symmetry 

Among the point groups, the reflection symmetries can 
be distinguished in two classes by the relative positions 
between mirror plane and rotation axis k z . In the first 
class, k z axis parallels the mirror plane, and the system 
doesn’t preserve inversion symmetry which correspond 
to point groups ^ 2 ^, CW, 64 ^, and T^. In these 


Mirror operator is an inversion operation followed by 
a C 2 rotation whose rotation axis perpendicular to the 
mirror plane. It should satisfy the following constraints: 
(1)[M, T] = 0, (2 )M+M = 1, (3 )MM = e**. If we 
set mirror plane as yz- plane or xy- plane, the reflec¬ 
tion symmetry operator will have the form: (A )M k = 
±To 0 i(?k, (B)Mfc = ±r z 0 icr k , (C )M k = ±r x (g) icr k and 
(D )M k = ±ir y 0 i(Jh{k = £, 2 , see Appendix), among 
which (C, D) represent glide mirror symmetry. For ex¬ 
ample, with the basis containing only equivalent sites [see 
FIG. (tj(ai)] |P/i+t), |-FU_t), \Pa-i), \Pa+i), a reflection 
operation along yz- plane writes M x = —ro 0 icr x . And 
if we set the basis as | Pa+i), I Pb+i), ~\ p A-t ), -j-Ps-t), 
where A/B are two inequivalent sites for the same kind 
of atom [FIG.[l](a 2 )], then the reflection operation along 
yz- plane with an additional interchanging of A and B 
sites writes M x = r x 0 ia x (see Appendix). Note that 
the basis in our framework is the single atom basis after 
considering spin-orbit coupling(SOC). 

Owing to the TRS and rotational symmetry, the 
Hamiltonian in k z axis has been constrained as Eq. 0- 
After we impose reflection symmetry, the H(k z ) can 
be further constrained and will reveal Dirac semimetal 
phase. As shown in Appendix, different mirror opera¬ 
tor will hold different type of Dirac semimetal. When 
the mirror operator is M k = tq/t z 0 a k (k = x,z\ only 
& 30 {k z ,m) survives. a^(k z ) is an even function respect 
to k zi so a 3 o(k z ,m) ~ Mq — M\k 2 z for the leading order. 
Tw o Dirac points will emerge in the k z axis at k z = 
± a/ Mq / Mi . This kind of Dirac semimetal is created 
through the two bands accidentally crossing each other 
when the conduction band and valence band have differ¬ 
ent rotation eigenvalues (p ^ q). It can be understood as 
a phase between normal insulator and weak topological 
insulator [42l 143] /topological crystalline insulator [44l [45] 
[FIG. |TJ (b)] - When mirror operator is M k = iT y /r x 0cr/c, 
only CL 3 s(k z ,m) survives which is an odd function respect 
to k z . The leading order is a 33 = M 2 k z . There is one 
Dirac point at TRIM generated by band crossing (TBC). 
Under this scenario, the Dirac points are created and sta¬ 
bilized by the crystalline symmetry [FIG. H(c)(d)]. We 
will use tables below to show the physical properties of 
all kinds of Dirac semimetals. 
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FIG. 1. (a) There are four kinds of reflection symmetries in space groups: the mirror plane parallel to rotation axis as shown 
in ai and a .2 and the mirror plane perpendicular to rotation axis as shown in a 3 and a.4. In both cases, for normal mirror 
plane ai and a 3 , there is only one equivalent site. As for glide mirror plane a .2 and a. 4 , which is a reflection symmetry with 
a t translation, there are two inequivalent sites in the lattice. (b),(d) correspond to the Dirac points of type ABC and TBC 
in TABLE [I] and TABLE [H] respectively. (b)The phase transition determined by a control parameter m and Dirac semimetal 
created by ABC. When m is in a proper range, two Dirac points show up on k z axis. This phase lays between two gapped 
phase such as normal insulator and weak topological insulator (WTI) and topological crystalline insulator (TCI). (c)The Dirac 
semimetal phase is protected by crystalline symmetry. In this case, Dirac point can be found at TRIM. (d)When the conduction 
band and valence band have the same rotation eigenvalue, i.e. p = q, they will never cross each other because of strong level 
repulsion. 


II.2. Classification table 

If we impose reflection symmetry, TRS and rotational 
symmetry to the Hamiltonian in k z axis, we will have the 
criteria for basis’s rotational eigenvalues to obtain Dirac 
semimetal phase. The classification of Dirac semimetals 
when the mirror plane is parallel to k z axis is shown in 
TABLE [T| 62 and 63 rotational systems can only gen¬ 
erate Dirac points via TBC. 64 and Cq symmetries can 
protect Dirac semimetal phase in the presence of all those 
four reflection symmetries. 

The classification of 3D Dirac semimetals when the 
mirror plane is perpendicular to k z axis is demonstrated 
in TABLE |II| Inversion symmetry can emerge through 
P = 62 where C 2 is a two-fold rotation along k z axis 
in systems preserving 62 , 64,66 symmetry. Yang and 


Nagaosa [38] have already considered unitary inversion 
operator P = tq,t x ,t z . We show that in TABLE |ll| the 
same results hold when inversion operators are unitary. 
However, the mirror operators can also generate the an¬ 
tiunitary inversion operator P = ir y . P = ir y is an 
inversion operator with a translation in nonsymmorphic 
space group. It will produce a phase factor after being 
applied twice. One special case is the D^h group, which 
does not involve inversion symmetry. As shown in TA¬ 
BLE^ only M z = ir y ® o z can protect a Dirac point 
at TRIM with the combination of rotation eigenvalues 
p = 2 , q = 0 or p = 2 , q = 1 . 
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TABLE I. Classification table when k z is parallel to mirror plane. Dirac semimetals can be obtained in the systems with reflection 
and C n rotational symmetry. Here we choose k z as the rotation axis and assume that yz- plane is a mirror plane, (p, q , r, s) can be 
regarded as a set of orbital angular momentum in z direction and j is the total angular momentum. For example, for (p, g, r, s) = 

• 7tT • 5 7T • 17T • 3tT . 1 7T *37T -17T ‘ 3 7T 

(3, 2, 0,1) in C 4 system, the C 4 rotation eigenvalues are (e 1 ^, e 1 ^, ), which equal For 

compact presentation, we assume n/2 < q < p < n and consider the equivalence between {p, r} and {(/, s}. The leading order 
of f±,g±,go + g z ,go — g z are shown in the table. Each term should by multiplied by an coefficient function of k z respecting 
to the parity of aij(k) when constructing the elements of the Hamiltonian. The bulk Dirac points are obtained through ABC 


when M x = ±to/t z 
group /4icd[U . 


) cr x and they are obtained through TBC when M x = ±T x /ir y 


The material CdaAs 2 belongs to space 


M x 

Cn 

(p,q,r,s) 

Total j 

f± 

g± go + gz 

go - g z 

Dirac 

Type 

Materials 

Dispersion in 
k y direction 

±to/t z 

C 2 

- 

- 

- 

- - 

- 

- 

- 

- 


c 3 

- 

- 

- 

- - 

- 

- 

- 

- 


c 4 

(3,2,0,1) 

( ± i> ± §) 

k+ 

k± k- 

+ 

ABC 

Cd3 As 2 (/4i cd) 

Linear Dirac 


c 6 

(5,4,0,1) 

(±K±§) 

k+ 

k 2 . k- 

ki 

ABC 


Linear Dirac 


c 6 

(5,3,0,2) 

(±K ± D 

k\ 

k 3 ± k- 

+ 

ABC 


Linear Dirac 


c 6 

(4,3,1,2) 

(±i±§) 

k+ 

k\ k 3 + 


ABC 


Linear Dirac 

£ 

-H 

c 2 

- 

- 

- 

- - 

- 

- 

- 

- 


c 3 

(2,1,0,1) 

(±K±i) 

k+ 

k+ k- 

k+k- 

TBC 


Linear Dirac 


c 4 

(3,2,0,1) 

(±K±§) 

k+ 

k 2 ± k- 

k+ 

TBC 


Linear Dirac 


c 6 

(5,4,0,1) 

(±K±D 

k+ 

kl k- 

ki 

TBC 


Linear Dirac 


Cq 

(5,3,0,2) 

(±§,±!) 

k 2 + 

k 3 ± k- 

+ 

TBC 


Linear Dirac 


c 6 

(4,3,1,2) 

(±i,±f) 

k+ 

k\ k± 

+ 

TBC 


Linear Dirac 

±ZTy 

C 2 

(U,0,0) 

(±KH) 

k+k- 

k+ k± 

k± 

TBC 


Linear Dirac 


c 3 

(2,2,0,0) 

(±K ± D 

k+k- 

k- k- 

k- 

TBC 


Linear Dirac 


c 3 

(2,1,0,1) 

(±i±§) 

k+ 

k+ k- 

k+k- 

TBC 


Linear Dirac 


c 4 

(3,3,0,0) 

(±§,±§) 

k+k- 

k- k- 

k- 

TBC 


Linear Dirac 


c 4 

(3,2,0,1) 

(±i±§) 

k+ 

k 2 ± k- 

k+ 

TBC 


Linear Dirac 


c 4 

(2,2,1,1) 

(±§,±§) 

k+k- 

k+ k+ 

+ 

TBC 


Linear Dirac 


Cq 

(5,5,0,0) 

(±KH) 

k+k- 

k- k- 

k- 

TBC 


Linear Dirac 


Cq 

(5,4,0,1) 

(±i±§) 

k+ 

kl k- 

ki 

TBC 


Linear Dirac 


Cq 

(5,3,0,2) 

(±K±§) 

4 

k 3 ± k- 

+ 

TBC 


Linear Dirac 


Cq 

(4,4,1,1) 

(±f,±l) 

k+k- 

k 3 ± k 3 ± 

ki 

TBC 


Quadratic Dirac 


Cq 

(4,3,1,2) 

(±§,±i) 

k+ 

k <2 + k± 

+ 

TBC 


Linear Dirac 


Cq 

(3,3,2,2) 

(±i±f) 

k+k- 

k+ k+ 


TBC 


Linear Dirac 


II.3. The coexistence of Dirac points and Weyl line 

nodes 

For now we can construct Dirac points in the sys¬ 
tems preserving reflection symmetry but without inver¬ 
sion symmetry. It is known that breaking inversion sym¬ 
metry can produce Weyl points or Weyl line nodes. So 
we expect in some cases Dirac points and Weyl line 
nodes can exist simultaneously with SOC. We illustrate 
this distinctive properties in materials with C± v point 
group like CdsAs 2 through k • p perturbation method 
(see Appendix). Note that we choose mirror operator as 
m x = r z 0 ia x with p=3,q=2,r=0,s=l corresponding to 
the 3rd row of TABLE |U The detailed calculations are 
shown in Appendix. By choosing some proper param¬ 
eters, the band structure shown in FIG. [2] (a) exhibits 
two Dirac points along T — Z direction and two Weyl 
line nodes in yz -plane and xz-plane. FIG. [2j (b) shows 
the phase transition between different topological phases. 


Obviously the phase transition of Weyl line nodes is in¬ 
dependent of the emergence of Dirac points in k z direc- 
tion[FIG. [2j (c)]. The simultaneous appearance of Dirac 
points and Weyl line nodes indicates a new class of topo¬ 
logical phase with time reversal symmetry and SOC. 

As illustrated in FIG. UK a) this k • p Hamiltonian can 
have different topological phases. These phases depend 
on the SOC term Dq and inversion breaking term Bq. 
FIG.[2j(b) shows the phase diagram respecting to Bq and 
D 0 . There are two different phases when changing pa¬ 
rameters. In the gray area, the systems break inversion 
symmetry but can not protect Weyl semimetal. In the 
dark blue area, the conduction band and valence band 
will cross to form a Weyl line node with Aq = —0.06eV\ 
When Ao = —0.00922eV the dark blue area extend to 
the light blue area . When the inversion breaking term 
B 0 « 0, the system will not have Weyl line nodes phase 
on the Dq axis, but it can protect Dirac semimetal phase 
just like in CdsAs2(/4i acd). This phase diagram is inde- 
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TABLE II. Classification table when k z is perpendicular to mirror plane. Dirac semimetals can be obtained in the systems 
with C n rotational symmetry and reflection symmetry. Here we choose k z as the rotation axis and assume that xy-plane is a 
mirror plane. (p,q,r,s) can be regarded as a set of orbital angular momentum in z. For example, for (p,q,r,s) = (2, 0,1, 3) 

• 5 7T • 17T • 37T • ll T • 3 IT • \ TV • 3 IT ’ 1 7T 

in Ca system, the Ca rotation eigenvalues are , e 1 ^), which equal j i s the total 

angular momentum. For compact presentation, we assume q < p < n and consider the equivalence between {p,r} and {q, s}. 
The leading order of f±,g±,go + g z , go — 9 z are shown in the table. Each term should be multiplied by an coefficient function 
of k z respecting to the parity of aij{k ) when constructing the elements of the Hamiltonian. The bulk Dirac points are obtained 


through ABC when M z = -p.ro/ t z 
belongs to space group /4i acd. 


and they are obtained through TBC when M z = ±r x /ir y 0 a z . The material Cd 3 As 2 


M z C 2 P 

Cn 

(p,q,r,s) 

Total j 

/± 

9± 

-H 

0 

Dirac 

type 

Materials 

Dispersion 
in k y 

To/r z To 

to/t z 

c 2 

- 

- 

- 

- 

- 

- 

- 

- 



C A 

(2,0,1,3) 

(±§>± 5 ) 

*± 

k+ 

0 

ABC 

CdsAs2(/4i acd) 

Linear 



c 6 

(2,0,3,5) 

(=•=§>=•=§) 

4 

ki 

0 

ABC 


Quadratic 



c 6 

(3,1,2,4) 

(±§>±§) 

4 

k- 

0 

ABC 


Linear 



c 6 

(4,0,1,5) 

(±|.±5) 

kl 

k- 

0 

ABC 


Linear 

r z 

Tz/To 

c 2 

- 

- 

- 

- 

- 

- 

- 

- 



C 4 

(1,0,2,3) 

(±f>±§) 

k+ 

k 2 ± 

0 

ABC 


Linear 



c 6 

(1,0,4,5) 

(±|.±5) 

k+ 

k 2 + 

0 

ABC 


Linear 



C 6 

(2,1,3,4) 

(±§,±§) 

k+ 

kl 

0 

ABC 


Linear 



c 6 

(3,0,2,5) 

(±§>± 5 ) 

kl 

kl 

0 

ABC 


Quadratic 

il~y To 

ITy 

c 2 /c 6 

- 

- 

- 

- 

- 

- 

- 

- 



Ca 

(2,0,1,3) 

(±§>±1) 

k 2 ± 

k- 

k + 

TBC 


Linear 

T z 

T x 

c 2 

(1,0,0,1) 


k± 0 

k± 

TBC 

Distorted spinels 46 

Linear 



Ca 

- 

- 

- 

- 

- 

- 

- 

- 



c 6 

(3,0,2,5) 

(±|i±|) 

ki 

0 

k + 

TBC 


Linear 



Ce 

(4,1,1,4) 

(±|,±l) 

ki 

0 

kl 

TBC 


Cubic 

T x To 

T x 

C 2 /C 6 

- 

- 

- 

- 

- 

- 

- 

- 



Ca 

(2,0,1,3) 

(±§.±§) 

k± 0 

k+ 

TBC 

Bi02 

Linear 

T z 

ITy 

C 2 /C 4 

- 

- 

- 

- 

- 

- 

- 

- 



c 6 

(3,0,2,5) 

(±§>±1) 

k 3 ± 

ki 

k + 

TBC 


Linear 


TABLE III. Classification table of topological phase for D^h- In the presence of C 3 symmetry, Dirac semimetals can only be 
obtained in the systems with M z = ir y 0 a z . H(k z ) is the Hamiltonian on k z axis after considering reflection symmetry but 
without constraints by rotation symmetry. 


Mirror 

operator 

H(k z ) 

Cn 

Constraint 

Possible 
(p, q, r, s) 

Total j 

Dirac 

type 

TO 

hqt — CL10T x + CL 2 3 Ty + CL 3 oT z 

= a 0 i - ia 0 2 + (ail - iai 2 )T x + (a 3 i - ia 3 2)T z 

p^q^r 

q^s 

— 

— 

— 

T x 

-f = CLioTx + 0-20 Ty + CL33T z 

hfi = aoi — iao2 + (an — ia\2)T x — i(a>2i — icL22)T y 

p^q^r 

Q ¥= s 

— 

— 


iTy 

h n = a 13 T x + a 2 3T y + a 3 3T z 
— aoi — iao 2 


(2,0,0,2) 
(2,1,0,1) 

(±§>±§) 

(±§,±i) 

TBC 

TBC 

T z 

h-tt — ai3T x + a 2 0Ty + Cl 30 Tz 

/in = aoi - iao2 - i(a 2 1 - ia 2 2)T y + (a 3 i - id 32 )T z 

v + q + r 

q^s 

— 

— 



pendent of the creation of Dirac semimetal phase in bulk 
band structure. 

Numerical calculation shows that the two crossing 
bands of Weyl line nodes in blue areas have different 
mirror eigenvalues, which protects the gap closing on the 
Weyl line nodes. 

Weyl line nodes are protected when the crossing bands 


have different mirror eigenvalues due to the absence of 
level repulsion. As is shown in FIG. |3j(c)(d), Weyl line 
nodes may emerge on various kinds of mirror plane. 
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FIG. 2 . (a) The electronic structure of Ca v system wit h SOC. The Dirac points and Weyl line nodes can exist simultaneously. 
The parameters are listed in the Appendix TABLE |vT| (b)Phase diagram with respect to Bo, Do and Aq of Ca v system results 
from the Hamiltonian [Equation ( |11| )] of k ■ p perturbation method. For all other parameters we fix them using the values fixed 
for Cd 3 As 2 [E](see Appendix TABLE [vT| . The inversion breaking term is Bo and SOC is described by Do . Weyl line nodes exist 
in the blue areas, among which the dark blue area represents the situation when the parameter Ao = 0 . 06 eV(Cd 3 As 2 )[l 8 ] and 
the light blue area (partially covered by the dark one) for Ao = 0.00922eV fHgTel [47j. (c)Schematic diagram of the distribution 
of Dirac points(yellow points) and Weyl line nodes(red or blue circles) in Brillouin zone. Dirac points locate on the interception 
of four Weyl line nodes. 


II.4. Dirac and Weyl line nodes in mirror plane 

Recently proposed by Fang, et al|23], Dirac line nodes 
in mirror plane emerge in systems with inversion sym¬ 
metry, time reversal symmetry and reflection symmetry. 
The combination of inversion operation and reflection 
operation generates C 2 rotation perpendicular to mirror 
plane, which is covered by our framework. Due to the 
combination of inversion symmetry and TRS, two bands 
related by TRS actually stick together and make up a 
two-fold degenerate band. 

When we study the Hamiltonian in the mirror plane, 
the mirror symmetry give MH{k x ,k y )\k z =oM~ x = 
H{k x ,k y )\k z =o- Therefore, if the conduction (valence) 
band consists of two bands with different mirror eigen¬ 
values, the bands with the same (positive or negative) 
mirror eigenvalue from the conduction and valence band 
will have strong level repulsion and open up a gap (FIG. 
[3j (a)(b)). Otherwise when TRS-partner bands have the 
same mirror eigenvalues, the conduction band will have 
a different mirror eigenvalue from the valence band, and 
the level repulsion will be relaxed. The band crossing 
between them will create Dirac line nodes in the mirror 
plane. Thus Dirac line nodes in mirror plane can only be 
protected by the nonsymmorphic M z = ±ir y (&i(j z (FIG. 


the crossing bands have different mirror eigenvalues. As 
is shown in FIG. |3j (c)(d) and TABLE IV, Weyl line 


U(a)(b) and TABLE |lV). 

Different from Weyl line nodes in HgCr2Se4[24| and 
(CdO )2 (EuO) 2 ,26j, we construct Weyl line nodes in the 
presence of TRS. Weyl line nodes are protected when 


nodes may emerge on various kinds of mirror plane be¬ 
cause of its two-fold degeneracy nature. Whereas Dirac 
line nodes in mirror plane can only be protected by 
M z = ±ir y (&i(j z which satisfies the requirement of mirror 
eigenvalues above. Note that DLN and WLN protected 
by mirror symmetry also apply to two-dimensional sys¬ 
tems, because in three-dimensional systems we just take 
a plane into consideration by fixing k z . 


III. CONCLUSION 

In this work we show that reflection symmetry can 
protect Dirac semimetal phase with or without inversion 
symmetry. We have classified Dirac semimetal in the 
systems preserving reflection symmetry, rotational sym¬ 


metry and TRS. TABLE IV can be referred to for an 
overall possible protection of semimetal phases by reflec¬ 
tion symmetry. There are two kinds of Dirac semimetals 
created via ABC and TBC. In C 2 and C 3 rotation invari¬ 
ant systems, Dirac semimetals can be created only via 
TBC. Whereas in C 4 and Cq systems, Dirac semimet¬ 
als can be created via not only ABC but also TBC. 
We also found that in C^v point group system, Dirac 
semimetal phase can coexist with Weyl line nodes. Fi¬ 
nally we show that DLN in mirror plane can be pro¬ 
tected only by M z = ir y ( 8 ) icr z and WLN can be pro- 
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/ 

\ 


Negative mirror 


eigenvalue 


Positive mirror 


eigenvalue 

\ 

J 




FIG. 3. Dirac or Weyl line nodes protected by reflection sym¬ 
metry. Only iT y 0 g z can protect Dirac line nodes while all 
kinds of mirror plane can protect Weyl line nodes. (a)(b)[23] 
are inversion preserving systems where Dirac line nodes may 
emerge and (c) (d) are inversion breaking systems where Weyl 
line nodes may emerge, (a) For mirror operators other than 
iT y ®icr z (e. g. t x ®icr z ), two bands with different mirror eigen¬ 
values compose conduction (valence) band. Apart from two 
Dirac points, a gap between conduction and valence bands will 
open due to level repulsion, (b) For mirror M z = ir y 0 icr z , 
conduction (valence) band consists of two bands with the 
same mirror eigenvalues. There isn’t level repulsion between 
conduction and valence bands and the crossing of two bands 
generates a Dirac line node. The existence and classification 
of DLN by different reflection symmetries are shown in TA¬ 
BLE |Tv| (c)(d) For various kinds of mirrors (e.g. ir y 0 ia x 
and r x 0icr x ), Weyl line nodes can be protected as long as the 
crossing bands have different mirror eigenvalues and so level 
repulsion doesn’t happen (TABLE |IV| ). 


Zihao Gao acknowledge financial support from Yat-sen 
school, Sun Yat-sen University. 

V. APPENDIX 

V.l. The constraint on Hamiltonian by TRS and 
rotational symmetry 

In this section we first constrain the Hamiltonian with 
TRS and rotational symmetry. Then we impose ro¬ 
tational symmetry C n to the whole system to deter¬ 
mine the leading order of the elements in h^(k) and 
The Hamiltonian and the basis we use here 
are the same as those in the main text. Time reversal 
symmetry will give constraints to the Hamitonian H{k ): 
H(—k ) = TiL(/c)T -1 , where the time reversal operator 
T - icr y K. We will have [38] 

H= ( h n$) h n (k) \ 

v h *n(-k J 

At the same time, we know the par¬ 
ity of each coefficient with respect to mo¬ 
mentum ^01,02,03,11,12,13,20,31,32,33 ( k') = 

~a01,02,03,ll,12,13,20,31,32,33(&) anc ^ a 00,10,21,22,23,30(~ k) = 
&oo,io, 21 , 22 , 23,30 (&)• We set the rotation axis as k z and 
choose the eigenstates of the rotation operator C n as the 
basis of matrices. Then the matrix representation of C n 
is 


TABLE IV. The possible protected semimetal phases by re¬ 
flection symmetry. The mirror operators with cr x stands for 
mirror parallel to rotation axis and cr z for mirror perpendic¬ 
ular to rotation axis. Dirac line nodes in mirror plane is de¬ 
noted by DLN. 


Mirror opertor 

ABC TBC WLN DLN 

T 0 0 icr x 

V 


V 

r x 0 icr x 


V 

V 

ir y 0 icr x 


V 

V 

r z 0 ia x 

V 


V 

To 0 icr z 

V 


V 

r x 0 icr z 


V 

V 

iT y 0 ia z 


V 

V V 

t z 0 icr z 

V 


V 


Cn — diag[dp, Qqi ^s] 


e ilf(P+l/2) 0 0 

0 e i ^(«+ 1 / 2 ) 0 

0 0 e^e+V 2 ) 

0 0 0 


27r( 


i+p+q | p—g 





0 \ 
0 

0 

e^wU+i/ 2 ) J 


( 2 ) 


where p, g, r, s E {0,1,..., n — 1} and can be regarded as 
orbital angular momentum of different states. In general, 
C n commute with TRS [C n ,T] = 0, thus p and r, q and 
s are related by: 


tected by any mirror operator. These new classes of Dirac 
semimetals in inversion breaking and preserving systems 
can guide the search for novel materials with exotic quan¬ 
tum properties [48]. 
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CM.p — Clr? ^ s 

exp [i—(p + r + 1)1 = l,exp [i — (q + s + 1)1 = 1 ^ 

n n 

Next we derive the constraint relations between rota¬ 
tional symmetry and elements of the Hamiltonian. The 
2x2 block Hamiltonian /qq, can be expanded in the 
following wav(49j: 

h n (k) = f 0 (k) + f+(k)T+ + /+(fc)r_ + f z (k)r z 
ftf-lA) = 9o(k) + g+(k)T+ + g-(k)r- + g z (k)r z 

























where r± = r x db ir y , fo :Z are real functions and 
f+,go,g z ,g± are complex functions. Then the rotational 
symmetry C n H(k±,k z )C~ 1 = H(k±e ±l27T / n , k z ) gives 
the constraints of elements of the Hamiltonian: 

fz{k±,k z ) = f z (,k ± e ±l2n/n ,k z ) 

f+(k±,k z ) = f + (k±e ±l2n/n ,k z ) 
go+z(k±, k z ) = g 0+z {k±e ±l27r/n , k z ) 
9o-z(k±, k z ) = g 0 - z (k ± e ±i2 ^ n , k z ) 
9±{k±, k z ) = g±(k±e ±l2 * /n , k z ) 


exp 

exp 

exp 

exp 


.271- 

i (p — q) 

n 

'.271r 

i — (p — r) 
n 

.271- ' 

i —(9 - s) 
n 

.27r 

i-(q-r) 

n 


(5) 


where fc± = A: x ± ik y: g 0 ± z = g 0 ±g z . 
On axis, these constraints become, 


exp 

exp 

exp 


’. 2 tt ' 

i—(p-q) 

n 

.271- . 

i — (p — r) 
n 


.271- 

7 — (q - s ) 
n 


exp 


• 2 tt ^ 

1 [q — r) 


fz(k z ) = f z (k z ) 

U(k z ) = u(k z ) 

9o+z(k z ) = 9o+z {kz) 
go-z(k z ) — 9o—z{kz) 
9±(k z ) = g±(k z ) 


( 6 ) 


After considering all other symmetric relations (e.g. re¬ 
flection symmetry in V.3), if non-diagonal /, g terms 
aren’t eliminated on k z axis, the corresponding (p,q,r,s) 
pairs (/+ to p,q, g 0 + g z to p, r , g 0 - g z to q,s, g± to q, r ) 
should be unequal to obtain Dirac points. For example, 
if /+ exists on k z axis, we should apply p ^ q or a gap 
will open on k z axis. 

In order to get the dispersion relation near Dirac 
points, we should change the /, g terms into a more ex¬ 
plicit form. The matrix elements can be expanded as 
polynomial j49y 


f(k+,k-) = A nin2 k^k n J. 


(7) 


ni,ri2 


Combined with the constraint relations ([ 5 |, we obtain 
e i2 ^ p -^/ n f(k+,k_) = f(k+e i2n / n ,k-e~ i27r/n ) 


= J2 ex P 


i—(ni - n 2 ) 
n 


A h ni h n<2 

/ 1 nin 2 A/ -|- l% — 


(8) 


where A niTl2 is an complex coefficient. To satisfy above 
equations, the phase factors must cancel each other, i.e. 
n\ — n 2 = (p — q) mod n. We choose the leading order 
terms to complete TABLE [I] and TABLE [II| For example, 
in C 4 system with p = 3 and q = 2, the constraint relation 
of /+ term is, 

6XP f + ( k±,kz ^ = f+^k± e±l ^ /n ,k z ) (9) 


Obviously the leading order term is given by n\ = 1, ri 2 = 
0 , and so we can replace the /+ term by Ai ? o&+, neglect¬ 
ing the higher order terms. 


V.2. Details about the classification of reflection 
symmetry operators 

We study the general classification of reflection sym¬ 
metry operators. The matrix representation of reflection 
symmetry can be decomposed to orbital space and angu¬ 
lar momentum space. In spin space, mirror operator is 
a two-fold rotation perpendicular to the mirror plane h: 
R n/2 (n) = e -* cr i- n7r / 2 5 where j indicates the half-integer 
spin momentum. At the same time, mirror operators 
should satisfy the following constraints: (1)[M, T] = 0, 
(2)M + M = 1, (3 )MM — e The reflection symme¬ 
try operators can take four possible forms: (A )M k = 
±r 0 (8) icr k , (B)M/c = ±r z (8) icr k , (C )M k = ±r x (8) io k and 
(D)Mfc = ±ir y <g> icr k (k = x,z, see Appendix), among 
which (C, D) represent glide mirror symmetry. 

As is shown in FIG. [I], physically there are two classes 
of mirrors, the mirror parallel to rotation axis and the 
mirror perpendicular to the rotation axis, respectively 
correspond to k = x and k = z above. 

We can now give the physical interpretation of these 
mirror operators below, where A/B represent two in¬ 
equivalent atom sites (see FIG. 0- \P±) = \P X )±i\Py), 
and t, i represent spin: 

1. We set the four basis as 

(l-FU+t); l-P-A-t)) \Pa~i), l-FU+ 4 .))- Then the 

yz- plane mirror operation interchange P + ,P_ and 
t, I, i.e. the mirror operation M x : | p A+\) 

— i\PA-t)i \ p A+t) — i\PA-i), \ p A-t) 

—i\P A +l)i | Pa-]) -> -i|PA+t)- The matrix repre¬ 
sentation of mirror operator is M x = —To 0 icr x . 

2. We set the four basis as 

(I Pa+i), Pb+i),-\Pa t)- ~\ p B-t))- A S lide P lane 
operator can have the following transformation 

M x : | Pa+i) -i\Pb -t)» I Pb+i) —i\PA- 1)» 

\PA-t) -> ~i\ p B+i),\PB-\) -> ~A p a+i )• There¬ 
fore the matrix representation of mirror operator 
is M x = =b r x 0 i(J x - 

In general basis for this mirror operator can be 
constructed as following: There are two inequiv¬ 
alent sites (A and B) and the distance between 
them is t. This transformation can be provided by 

defining: |P A ±ta)) = e±i ^ Ku A±m)A P B±Ui)) = 
e ±iMr-K UB ± where K denote the point in Bril- 

louin zone, ua±(Mt + ^)e ±lt ' K = ub±(v) and 
UA{r+Mt+t) = ua{t)> M z : r Mr+t , where M 
is a symmorphic mirror operation acting on k space 
and K • (Mt -f t) — 27m. For Bi02 of space group 
7^227, the Dirac point appears in the X point of its 
Brillouin zone with reflection symmetry belong to 
this case. 
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3. We set the four basis as 

( \PA+t),\ p A+i),\ p A-i),-\PA-t )) yz- plane mirror 
operation M x : \Pa+i) -i\Pa-\), |TU+t) 
—i\PA-i), l-FU-t) -i|PW,|-FVt> ~^ 

—i|PA + ^) has a matrix form M x = —r z <S> icr x . 

4. We set the four basis as 

(I Pa+i), I Pb+i), ~\Pa-i),~\P B - t))- The 

mirror operation on xy- plane writes M z : 

I Pa+i) i\Pb+i), | Pb+i) —^ -£|Pa+^), I^A-t) 

—i|P#_^), |Pg_^) —>• i\Pa-^): with a matrix 

representation M z = ir y ® icr z . There are two in¬ 
equivalent sites(A and B) and the distance between 
them is t. This transformation can be provided by 
defining: |-P4±t(t)) = e ± "' K u A ±^i), |-Pb±t 4)> = 
e ±iMr-K UB ± where K denote the point in 

Brillouin zone, u A ±(Mr + t)e ±lt ' K = ub±(t) and 
1^4 (r + Mt + i) = u A (r). M z : r Mr + £, where 
M is a symmorphic mirror operation acting on k 
space and K • (Mt + t) = (2n + l)7r. 

V.3. The constraint on Hamiltonian by reflection 
symmetry and rotational symmetry 

Now we show how the Hamiltonian in k z axis is con¬ 
strained by reflection symmetry and rotational symme¬ 
try: 

Mirror plane parallel to k z axis: First we set the mir¬ 
ror plane as yz plane and the mirror operator as M x . 
The combination of M x and n-fold rotation operation 
generate n-1 more mirror planes and therefore all these 
mirror operators can be denoted as M & = C^M x (k = 
0,1,..., n — 1). These mirror planes cross on k z axis 
and confine the Hamiltonian in k z axis. For each M&, 
there’s commutation relation M},H(k z )M^ 1 = H(k z ). 
For the basis satisfying constraint relations (|6|, when 
M x = to/t z ® ia x , the Hamiltonian becomes H(k z ) = 
diag[a 0 + a 30 ,a 0 - a 30 ,a 0 + a 30 ,a 0 - a 30 ]; when M x = 
TxliT y ®iG x , the Hamiltonian becomes H(k z ) = diag[ao + 
a 33 ,ao — a 33 ,ao + a 33 ,ao — a 33 ]. It is notable that at 
TRIM all a terms with odd parity eliminated due to the 
relation [ H , T] = 0 and only aoo, nice a 3Ch a 2 i, &22, a23 sur¬ 
vive. Therefore, the constraints on (p, g, r, 5) are lessen in 
TBC. All the results of the physical properties of Dirac 
semimetal are shown in TABLE [H 

Mirror plane perpendicular to k z axis: First we set 
the mirror plane as xy plane and the mirror operator 
as M z . Along the k z axis, the reflection symmetry is 
M z H(k z )M~ l = H(-k z ). Combined with C 2 ,C 4 ,C 6 
rotational symmetry, reflection symmetry can create in¬ 
version symmetry P = The rotation operator 

for orbital space is C n = i exp(y7r 1+ ^ +g ) exp(Pr^^r^) = 
i exp(y7r 1+ ^ +g ) [cos ^(p — q) + ir z sin ^(p — g)]. Two-fold 
rotation is C 2 = cos Oto + sin where 0 = § (p — g). 
If p — g = 0,2,4, then C 2 = ±to. These conditions 
can be satisfied when p = g for C 2 symmetry, when 
p = g or p = g + 2 for C4 symmetry, when p = g, 


p = g + 2 or p = g + 4 for (^6 symmetry. If p — g = 1, 3, 5, 
then C 2 = ±Tz. These conditions can be satisfied when 
p = 1, g = 0 for C 2 symmetry, when p = g + 1 or 
p = g + 3 for C4 symmetry, when p = g + l,p = g + 3or 
p = g + 5 for C6 symmetry. The symmetry constraints 
with TRS along k z axis are M z H(k z )M ~ 1 = TH(k z )T ~ 1 
and PiL(^)P- 1 = TH(k z )T~*. After [P,T] = 0 at 
TRIM and constraint relations |6| are applied, all sit¬ 
uations protecting Dirac points are shown in TABLE 
[H| holding the same result with Yang and Nagaosa|3U 
when P = zLtq , ±Ts, Try. Note that there are additional 
conditions for inversion symmetry PT^P -1 = ±r M (/i = 
03, 30, 33) |2H], the case for mixed mirror operator is thus 
eliminated. When P = ±ir y the system can also gen¬ 
erate Dirac points. For C 3 symmetry, the system does 
not have inversion symmetry. We only have symmetry 
constraints m z H(k z )m ~ 1 = TH(k z )T~ 1 and [P, T] = 0 
at TRIM and the result is shown in TABLE IIIII . 


V.4. k • p model 

In this section, we use the four band model to describe 
the effective Hamiltonian of material with point group 
C± v such as Cd 3 As 2 . We write the 4x4 effective Hamil¬ 
tonian generally as: 

3 

P e ff = ^ dij(k)Tij (10) 

i,j =0 

where Yij = r^aj. The basis of effective Hamiltonian are 
four angular momentum eigenstates | + ^),| + |),| — |) 
and | - |). 

Then we investigate the representations of T matrices 
and k polynomials with three symmetric operations: the 
four-fold rotation along Z axis C4, the vertical reflection 
m v and the dihedral reflection m& = C^m v . 

The matrix operation of symmetric operators are: 
(1 )Z7(C 4 ) = Pi(C 4 ) ® P|(C 4 ), (2 )U(m v ) = Ri(m v ) ® 

P| (rn v ), and (3 )U(rrid) = U ( m v )U (C 4 ), where Rj(C±) = 
ex.p(i^ja z ),Rj(m v ) = exp(i|jcr x ). 

The operators act on the k polynomials dij(k) as : 
(l)C7 4 . k x y /cy, k y y k x , k z y k z , (2 ^rfi v . k x y 
k x , k y —>• k y ,k z —>> k z , (3)rnd : k x —y k y ,k y —>> k x ,k z —>• 
k z . 

The representations of T matrices and k polynomials 
d(fc) are shown in TABLE [V| By assembling the T ma¬ 
trices and d(fc) with the same representation and time- 
reversal eigenvalue we obtain our effective Hamiltonian: 


FFeff(fc) = e(k)+ 



! A(k) -D(k)k+ 

B+k- 

—C(k) \ 

—D*(k)k- - A(k ) 

C(k) 

B_k+ 


m 

D(k)k- 

V -C*(k) B_k_ 

D*{k)k+ 

-A(k) / 


where e(k) = Eo~\-Eik + k--\-E2k^, A(k) = + 
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TABLE V. The representations of T matrices and k polyno¬ 
mials 


Reps 

T matrices 

d(k) 

T 

fi 

Too, T30 

1 ,kl + kl,k 2 z 

+ 

fi 

- 

k z 

- 

f 2 

r 0 3,r 3 3 

- 

- 

f 3 

r 22 

kl - kl 

+ 

f 3 

r 12 

- 

- 

f 4 

r 2 i 

kx ky 

+ 

f 4 

Tn 

~ 

- 

f 5 

(Tio, r 2 3) 

(k x k z , k y k z ) 

+ 

f 5 

(r 32 , Toi), (To 2 , T3i), (r 2 o, — Ti3) 

(k x , ky) 

- 


TABLE VI. Parameters for the 4x4 effective Hamiltonian. 
Most parameters can be referred to from Jeon, Sangjun, et 

al[E]. 


Eo(eV) 

-0.219 

Bo(eV A) 

5 

Ei(eV A 2 ) 

-30 

Bi(eV A) 

0 

E 2 (eVA 2 ) 

-16 

C 0 (eVA 2 ) 

0 

Ao(eV) 

-0.060 Ci(eVA 2 ) 

0 

Ai(eV A 2 ) 

18 

D 0 (eVA) 

-2.75 

A 2 (eV 2 A 2 ) 

96 

D\{eV A 2 ) 

0 

A 20 (eV 2 ) 

0.050 




y/A 2 k* + A |o, B ± = ±B 0 + B u C{k) = ^{k\ + k 2 _) + 
\C\(k\ — k'L), D{k) = Do + iD\k z , and k± = k x ± ik y . 

The Hamiltonian shows the possibility of coex¬ 
istence of 3D Dirac points and Weyl line nodes. By 
choosing proper parameters shown in TABLE |VI| (Only 
B 0 is different from CdsAs2) we can verify our state¬ 
ment through the band structure FIG. [2j (a). FIG. 
[21(a) is accomplished with the substitutions: ki -A 
sin (kiLi), k% -A A[l —cos (kiLi)\ for a periodic lattice, 

where L x = L y = a = 12.67A and L z = c = 25.48AjT7]. 

In order to solve Weyl line nodes, it is convenient to 
neglect the symmetric terms e(k) and some parameters 
B \, Co , C \, Di as is shown in TABLE [VT| Therefore, ex¬ 
plicit calculation gives the equation of Weyl line nodes: 


(B 2 0 - D 2 )k 2 x 



( 12 ) 


This is schematically shown in FIG.[2](c). 

In FIG.[2j the Co, Ci, D\ terms are neglected. B 0 
and Dq are two main parameters determining whether 
there is a gap or Weyl line nodes in BZ. The phase dia¬ 
gram of Bo and Do is shown in FIG. [2] (b). The critical 
line where Weyl line nodes appear and disappear is close 
to the line Do = Do, but Bo is actually slightly larger 
than Dq, which can be shown more clearly with a smaller 
Ao [see FIG. ||(b) the dark blue area]. On the critical 
line, the two bands pull apart and the Weyl line nodes 
eventually annihilate accompanied by opening a gap. 
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